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Spin-orbit coupled jes = 1/2 iridium moments on the geometrically frustrated fee lattice 


A. M. Cook^, S. Matern^, C. Hickey^, A. A. AczeP, and A. Paramekanti^’^ 
^Department of Physics, University of Toronto, Toronto, Ontario, Canada MSS 1A7 
^Institute for Theoretical Physics, Cologne University, 50937 Cologne, Germany 
® Quantum Condensed Matter Division, Oak Ridge National Lab, Oak Ridge, TN, 37831, USA and 
^Canadian Institute for Advanced Research, Toronto, Ontario, MSG 1Z8, Canada 

Motivated by experiments on the double perovskites La2ZnIr06 and La2MgIr06, we study the magnetism of 
spin-orbit coupled y'eff = 1/2 iridium moments on the three-dimensional, geometrically fmstrated, face-centered 
cubic lattice. The symmetry-allowed nearest-neighbor interaction includes Heisenberg, Kitaev, and symmetric 
off-diagonal exchange. A Luttinger-Tisza analysis shows a rich variety of orders, including collinear A-type 
antiferromagnetism, stripe order with moments along the {lll}-direction, and incommensurate non-coplanar 
spirals, and we use Monte Carlo simulations to determine their magnetic ordering temperatures. We argue that 
existing thermodynamic data on these iridates underscores the presence of a dominant Kitaev exchange, and 
also suggest a resolution to the puzzle of why La 2 ZnIr 06 , but not La 2 MgIr 06 , exhibits ‘weak’ ferromagnetism. 


Introduction. — Heavy atoms with strong spin-orbit cou¬ 
pling (SOC) and electronic correlations are predicted to form 
exotic quantum phases fT). Rare-earth ions with strong SOC 
on the frustrated pyrochlore lattice can yield local moments 
with unusual exchange couplings, leading to ‘quantum spin 
ice’, as in Yb2Ti207 IMI. Another exciting proposal is to 
realize the Kitaev Hamiltonian, with a spin liquid ground state 
and Majorana fermion excitations a, in iridium oxides with 
edge-sharing octahedra, such as the two-dimensional (2D) 
honeycomb iridates Na2lr03 and Li2lr03 EllH. Doping such 
Mott insulators has been predicted to lead to topological su¬ 
perconductivity lHOHIl- Experimentally, in both Na2lr03 
and Li2lr03, the spin liquid state is preempted by magnetic or¬ 
der HEM induced by interactions beyond the Kitaev model. 
Nevertheless, extensive work on these materials CMa, and 
3D harmonic honeycomb iridates /3,7-Li2lr03 ||2TI4301 . as¬ 
cribes their complex order to large Kitaev couplings. Kitaev 
interactions in the triangular iridate Ba3lrTi209 may lead to 
vortex crystals or gauge-like degeneracies EiHia. 

In light of these studies, we explore the following impor¬ 
tant issues. What kinds of phases does the Kitaev interaction 
support in 3D lattices with geometric frustration? Do experi¬ 
ments suggest dominant Kitaev interactions in any geometri¬ 
cally frustrated materials? Here, we address these questions in 
the context of ordered double perovskite (DP) compounds, a 
large class of materials with the chemical formula A2BB’06, 
where B and B’ ions occupy the two sublattices of a 3D cubic 
crystal. Metallic DPs such as Sr2FeMo06 ll34l are of great 
interest as half-metallic ferromagnets llT5l - l39l . Recent work 
on metallic DPs has examined the role of SOC on bulk spin 
dynamics ll40ll . and Chern bands in ultrathin films EUEl. 
On the other hand, DPs where B is an inert filled-shell ion, 
and B’ is a heavy Ad/5d ion, form Mott insulators with lo¬ 
cal moments on the frustrated fee lattice of B’ ions 
Our work is motivated by the recent synthesis of La2ZnIr06 
and La2MgIr06 ll55l . Structurally, both materials have nearly 
undistorted oxygen octahedra. A nominal valence Ir'*+ (5c?®), 
together with the strong SOC and larger spacing between Ir 
ions compared to perovskites, suggests that these materials 
behave as effective jeff=l/2 Mott insulators ll55l . 


In this Rapid Communication, we focus on the broad as¬ 
pects of magnetism in an ideal fee lattice, highlighting the 
rich physics of strong SOC in a canonical frustrated 3D lat¬ 
tice. Our key results are the following, (i) We show that 
even the nearest-neighbor symmetry-allowed Hamiltonian on 
the fee lattice, which includes Heisenberg, Kitaev, and sym¬ 
metric off-diagonal exchange couplings, leads to rich mag¬ 
netic phases such as collinear antiferromagnetism, stripes, or 
multimode spirals. Indeed, previous work 15^ has suggested 
that strong Kitaev interactions should be present in a large 
class of 2D and 3D lattices, including the fee lattice, but did 
not study the most general symmetry-allowed Hamiltonian, 
(ii) We find that strong SOC can also stabilize a regime of 
robust A-type antiferromagnetism (AFM), also called Type-I 
AFM, which is observed in neutron diffraction on Fa2ZnIr06 
and Fa2MgIr06 ll55l . Our results challenge the conventional 
wisdom which ascribes robust A-type antiferromagnetism in 
many fee magnets to further neighbor Heisenberg exchange 
lEIlIsll, and suggests that anisotropy due to SOC may be 
crucial in 5d oxides. Indeed, a recent ab initio study of 
Sr2CrSb06 1591 finds next-neighbor interactions are negligi¬ 
ble, < 5% of the first neighbor interactions, (iii) In certain 
regimes with A-type AFM, we uncover a residual accidental 
XY degeneracy of collinear states. Thermal order by disorder 
pins the moments along the Ir-O bond directions, (iv) We ar¬ 
gue that thermodynamic data on Fa2ZnIr06 and Fa2MgIr06 
lEa, i.e., their ordering pattern and small frustration param¬ 
eter, indicate a dominant antiferromagnetic Kitaev coupling. 
Microscopically, this may arise from the near-cancellation 
of Heisenberg interactions, from multiple Ir-O-O-Ir superex¬ 
change paths minil^lShl, and the smaller direct exchange for 
well-separated Ir atoms in the DP structure. We argue that a 
subtle difference in magnetic orders can reconcile ‘weak’ fer¬ 
romagnetism in Fa2ZnIr06 with its absence in Fa2MgIr06. 
These compounds thus realize a new class of ‘Kitaev materi¬ 
als’. Ultrathin films of Fa2BIr06, grown along {111}, could 
realize the triangular lattice AFM Kitaev model. 

Model. — To construct a minimal model on the fee lattice 
of Ir moments, we consider the ideal cubic DP structure, and 
focus on nearest neighbor terms which are expected to dom- 
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FIG. 2: Real space spin configurations in the layered A-type antifer¬ 
romagnetic states AFM A-I and AFM A-II. 


FIG. I: Phase diagram of the nearest-neighbor spin Hamiltonian as 
a function of Jk/Jh and F/ Jh, obtained using the Luttinger-Tisza 
(LT) method. The AFM states are A-type antiferromagnets, having 
ferromagnetic planes stacked antiferromagnetically along the third 
direction, with spins either pointing perpendicular to the FM plane 
(A-I AFM) or lying in the FM plane (A-II AFM). Stripe order at 
( 71 / 2 ,7r/2,7r/2) features moments pointing in the {111} or symme¬ 
try related directions. IC-1, IC-2, and IC-3 are incommensurate non- 
coplanar spirals; beyond the LT analysis, they are multimode states. 

inate. We appeal to symmetry arguments to write down all 
possible terms, based on the fact that the effective jeff = 1/2 
angular momentum operator is a pseudovector (axial vec¬ 
tor). Requiring invariance of the Hamiltonian under lattice 
rotational and mirror symmetries II 6 OI constrains the Hamil¬ 
tonian coupling nearest-neighbor Ir sites to be of the form 
H = iTn + Hk + Hob, 

Hh = Jh^2 

(rr'l 

Hk = Jk{Y.S^S^, + Y,S^S^, + ( 2 ) 

{rr’)^y (rr')a^ (rr')x^ 

Hob = rY, 

r 

- 5'^5'r-x-r!/) + {x,y^ y, z) + {x,y^ x, z)] (3) 

Here, (rr') denotes all first-neighbor pairs, while {rr')xy de¬ 
notes first-neighbors restricted to the ccy-plane (similarly for 
yz, xz). Hk is the Heisenberg term, Hk is the Kitaev in¬ 
teraction, and Hob is a symmetric off-diagonal exchange 
term. Antisymmetric Dzyaloshinskii-Moriya interactions are 
forbidden here by inversion symmetry. A dominant Jh < 0 
leads to ferromagnetism; this is incompatible with the order¬ 
ing observed in La 2 BIr 06 (B=Mg,Zn), so we assume Jh > 0. 

Luttinger-Tisza analysis. — To determine the pre¬ 
ferred magnetic orders, we use the Luttinger-Tisza (LT) 
method which considers the spins to be classical moments, 
and replaces the constant length spin vectors by uncon¬ 
strained vector fields (j)r. The classical spin Hamiltonian 
written in momentum space then takes the form Hih = 


Ek with 

/Ak + aCf \ 

M(k) = Ak + . (4) 

V - 75 - -75E Ak + aC^y 

Here, Ak = (cos kx cos ky + cos k^ cos kz + cos ky cos kz), 
= cos ki cos kj, and = sin ki sin kj, and we have 
defined a = Jk/Jh and 7 = T/Jh- Here, ki (with i = 
x,y,z) denote components of the momentum along the cubic 
Ir-O axes, and we have set the Ir-O-B bond length (B=Zn,Mg) 
to unity. Diagonalizing Hhc for Jh > 0 , and looking for 
the lowest energy eigenvalue in k, we find the rich variety of 
magnetic orders shown in Fig.[T] 

Magnetic orders. — The LT analysis yields collinear as 
well as spiral antiferromagnetic (AFM) states. We describe 
these phases below, and compare their energy with numerical 
simulated annealing results. 

A-I AFM: This is an A-type collinear AFM (also referred to 
as a Type-I AFM in the literature) which consists of ferromag- 
netically ordered spins in the cubic ab-plane layered antifer¬ 
romagnetically along the c-axis . The spins point along the 
c-axis, perpendicular to the ferromagnetic planes as shown 
in Fig. 1^ There are six symmetry related A-I AFM ground 
states, associated with a three-fold choice of the layering di¬ 
rection and a two-fold choice of the Ising AFM order. Al¬ 
though these are the lowest energy collinear states, there is 
an accidental classical degeneracy, where one can form mul¬ 
timode states leading to coplanar or even noncoplanar states 
with the same classical ground state energy. This degeneracy 
is expected to be broken in favor of collinear states by fluc¬ 
tuation effects, and our simulated annealing finds the above 
collinear states to be stabilized by thermal order by disorder. 

A-II AFM: This is also an A-type collinear antiferromagnet; 
however spins lie in the ferromagnetic planes as in Fig. In 
addition to collinear states, there are again multimode copla¬ 
nar or noncoplanar states with the same classical ground state 
energy; we expect and observe numerically that thermal fluc¬ 
tuations favor the collinear orders. However, the ground state 
energy is independent of the precise angle in the plane so that 
there is an accidental XY degeneracy of collinear states. Our 
simulated annealing results show that this degeneracy is also 
broken by thermal fluctuations, with ‘order by disorder’ fa¬ 
voring spins along the Ir-O bond direction. There are twelve 
symmetry related A-II ground states favored by fluctuations. 
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FIG. 3: Real space spin configurations in the collinear stripe state, 
showing moments pointing along the diagonal {111} and {Ill} di¬ 
rections for = (7r/2, 7r/2, 7r/2). 


arising from a three-fold choice of the layering direction and 
a four-fold choice of the spin axis. Remarkably, the A-IIAFM 
order persists even in the pure Kitaev limit with Jk > 0. 

Stripe: The collinear stripe state has spins pointing along 
the {111} and {Ill} directions arranged as shown in Fig. 
for (fca;, fcy, = ±(7r/2, 7 r/ 2 , 7 r/ 2 ); symmetry related orders 
are degenerate. Ordering with this wavevector is also referred 
to as a Type-II AFM. The ordering wavevector determines the 
direction of the spins, so that flipping one of the momentum 
components also flips the corresponding spin component; or¬ 
dering at ±(7r/2, — 7 r/ 2 , tt/2) leads to spins along {ill} and 
{III}. This leads to a total of eight ground states. 

Incommensurate Spiral (IC-1, IC-2): In these regimes, the 
LT analysis suggests an incommensurate coplanar spiral or¬ 
der with wavevector {kx,ky,kz) = (tt, Q, Q), and symme¬ 
try related equivalents. With a = Jk/Jh and 7 = T/Jh, 
minimizing the LT energy leads to Q = the 

transition into the AFM A-I state {Q — 0) happens when 
a = 2 | 7 |. However, we And that if we assume single mode 
ordering, the spins constructed in the IC-1 and IC-2 phases 
from the LT eigenvectors do not satisfy the constraint of con¬ 
stant magnitude. Our simulated annealing numerics show that 
the ground states in this regime are noncoplanar multimode 
spirals formed by superposing all six equivalent wavevectors 
(tt, Q, ±Q), (Q, tt, ±Q), (Q, ±Q, tt). 

Incommensurate Spiral (IC-3): In this regime, the LT ap¬ 
proach again suggests an incommensurate coplanar spiral or¬ 
der; however, the wavevector is of the form {kx,ky,kz) = 
{P, Q, Q). We have not found a simple closed form expres¬ 
sion for P, Q; however, they are obtained by minimizing the 
LT eigenvalue 

A = ( 4 - 1 - 0 !) cosPcos(3-l-(2-|-a) cos^Q — jsin^Q—VD (5) 

D = [o (cosP—cosQ) cosQ —7 sin^ Q] ^-F 87 ^ sin^ Psin^ <3 ( 6 ) 

Again, a single mode spiral does not satisfy the spin con¬ 
straint, and our simulated annealing numerics show noncopla¬ 
nar multimode spiral order in this regime. 

Monte Carlo results. — To complement the LT analy¬ 
sis, we have used simulated annealing numerics, which pre¬ 
serves the spin constraint, to And the classical ground states. 
Fig.Qa) compares the numerically computed ground state en¬ 
ergy per spin to the Luttinger-Tisza result, for Jk / Jh = 1.5 



FIG. 4: (a): Comparison of the ground state energy per spin Egs ob¬ 
tained within LT method (solid line) and simulated annealing (dots). 
(b),(c): Magnetic transition temperature Tc of the classical model 
(in units of JhS^, for spin length S) vs. T/Jh, obtained using 
Monte Carlo simulations for cuts through the phase diagram (Fig.[TJ 
at Jk/Jh = —1.0,-1-1.5. (d) Plot of the “frustration parameter”, 
the ratio of the TF = Tc}! -F I/S') to Qcw', the rescaling of Tc by 
(1 + 1/S) accounts for the classical S^ being replaced by the quan¬ 
tum S(S -F 1). The dark square shows the result at Jk/Jh —^ 00 . 


and varying F / Jh- The agreement between the two is excel¬ 
lent in the A-II AFM and Stripe states, where the collinear or¬ 
der is precisely recovered. Our result that the A-II AFM state 
appears even for large Jk differs from an earlier study ||5^ 
which proposed a spiral ground state based on a Luttinger- 
Tisza analysis which did not take into account thermal fluctu¬ 
ations and order-by-disorder. For IC-l/IC-2, the simulations 
indicate multimode order, and lead to an energy per spin (for 
36^ lattice) which is only slightly higher by < 2%. 

In order to determine the magnetic ordering temperature in 
the various phases, we have carried out Monte Carlo simula¬ 
tions on system sizes with up to 24^ spins. Fig. shows the 
magnetic Tc as determined from the specific heat singularity, 
along various cuts through the Luttinger-Tisza phase diagram. 
The Heisenberg limit in the absence of SOC {Jk = 0, F = 0) 
is the most fragile state with the lowest Tc ~ 0.44J//S'^; 
our results here agree with previous work on the fee Heisen¬ 
berg model ED, where thermal order by disorder leads to a 
nonzero Tc- The A-I AFM, A-II AFM, and stripe phases ap¬ 
pear most robust with high Tc, since SOC enhances the pin¬ 
ning of the moment direction. Thus, although the exchange 
interactions induced by SOC are frustrated on the fee lattice, 
the SOC nevertheless enhances Tc by favoring certain spin 
orientations, thus reducing the effects of thermal disordering. 

Comparison with experiments. — La2Znh06 and 
La2MgIr06 are A-type AFMs. Combined ab initio and neu¬ 
tron diffraction studies ll55l suggest that the Ir spins lie pre¬ 
dominantly in the ferromagnetic planes, viz. the A-II AFM 
state. This is consistent with Jk > 0 and |r| < Jk/2. Order 
by disorder pins moments along the Ir-O bond directions. 
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FIG. 5: Conjectured alignment of staggered octahedral rotations and 
the ferromagnetic planes in the A-II AFM state for La2ZnIr06 (top) 
and La 2 MgIr 06 (hottom), with the spins shown on Ir octahedra (yel¬ 
low). We have picked the a-axis as the direction along which the oc¬ 
tahedral rotations are staggered for the Ir octahedra, and shown only 
2 : = 0,1 planes. For La2ZnIr06, the FM planes are the a:t/-planes 
stacked antiferromagnetic ally along z, leading to a net ‘weak’ fer¬ 
romagnetic moment along —y, while for La2MgIr06 the FM planes 
are a: 2 -planes stacked antiferromagnetically along y leading to no net 
ferromagnetic moment. The uniform Ir octahedral tilts are unimpor¬ 
tant for this discussion and is not shown. 


The Curie-Weiss temperature of jeff = 1/2 moments on 
the ideal fee lattiee is Qcw = —{^Jh + Jk), independent of 
r. However, both La2ZnIr06 and La2Mgh'06 have a mon- 
oelinie P2i/n strueture, arising from small IrOe oetahedral 
rotations — an oetahedral rotation (p about the eubie c-axis 
whieh is staggered between adjaeent ab layers, and a global 
tilt about the eubie {110} axis. In the strong SOC limit, the Ir 
moments traek the oetahedral rotation, as shown for Sr 2 lr 04 
||8]|63- a high temperature expansion yields a powder aver¬ 
aged Ocw = —Jh— |(2J4/ +Jif)(l + 2cos2(/). If the axis 
along whieh the ferromagnetie planes are staeked in staggered 
fashion eoineides with the axis of the staggered oetahedral ro¬ 
tations, it leads to a net ferromagnetie moment « m sin p in 
the A-II AFM state, where m is the ordered moment. Equiv¬ 
alently, we may start with the Hamiltonian in the ideal eu¬ 
bie limit, and eonstruet the Hamiltonian for the ease with 
oetahedral rotations by making loeal unitary rotations on the 
j = 1/2 spins whieh induees Dzyaloshinskii-Moriya interae- 
tions, leading to an AFM with ‘weak’ ferromagnetism llsl l^ . 

In La2MgIr06, ab initio studies prediet a ‘weak’ ferromag¬ 
netie moment « O.S/is in the monoelinie P2i/n strueture; 
however, experiments do not deteet any fetTomagnetie mo¬ 
ment in the ordered phase. To understand this diserepaney we 
propose that the axis of the staggered oetahedral rotations and 
the staeking direetion of the fetTomagnetie planes are along 
orthogonal eubie axes (see Fig. |^, and ab initio results may 
have missed the eorreet ordering due to subtle energy differ- 
enees. This ean be tested if additional magnetie Bragg peaks 


ean be resolved using high resolution X-ray diffraetion. If we 
ignore SOC ( Jk = 0, T = 0), and note that « 9° from the 
struetural data is small, the measured Qcw ~ —24iT yields 
Jh ~ 8 iT. Our Monte Carlo simulations at Jh = 0, T = 0 
show Tc « OAAJhS^, consistent with previous work on the 
fee Heisenberg model 1^ . Heuristically replacing the classi¬ 
cal 5^ by S'(S' + 1) for quantum spins leads to a renormalized 
T/ = Tc(l -F S). This is a good approximation for the 3D 
cubic lattice S' = 1/2 Heisenberg model || 6 ^ . Here, on the 
fee lattice, with Jh = SK and S = 1/2, we find « 2.QK, 
much smaller than = 12A'. With T 7 F 0, but keeping 
Jk = 0, Tc hardly changes or even gets suppressed. This hints 
at a significant Jk > 0. Indeed, the “frustration parameter” 
/ = —Qcw/Tc, plotted in Fig. |^d) for T = 0, shows that 
recovering the experimentally observed small / « 2 needs a 
large Kitaev exchange Jk / Jh ^ 1- 

Thus, we suggest that a model with a dominant Kitaev term 
Jk > 0, perturbed by a weak Heisenberg exchange coupling 
Jh ^ Jk, is a good starting point to understand jeff = 1/2 
magnetism in La 2 MgIr 06 ; we estimate this dominant cou¬ 
pling Jk ~ 2AK. These estimates do not shed much light 
on the off-diagonal symmetric exchange since the powder av¬ 
eraged Qcw is independent of T, and Tc is not very sensitive 
to r (see Fig.j^c)). However |r| > Jk/2. is precluded by the 
observed order. Traditionally, in fee magnets, robust A-type 
order is ascribed to second-neighbor Heisenberg interactions 
Ell EH. For heavy oxides, however, our results show that the 
A-type AFM, and the small frustration parameter, is due to 
SOC-induced Kitaev interactions. 

In La2ZnIr06, there is a measured ‘weak’ ferromagnetic 
moment « 0.22/4^; thus, the axis along which the ferromag¬ 
netic planes are stacked in staggered fashion must coincide 
with the axis of the staggered octahedral rotations (see Fig.|^. 
Setting (p « 11°, consistent with structural data, we expect 
a moment « 0.19/iB, close to the measured value. This is 
smaller than the ab initio prediction « O.h/is. Based on the 
smaller T^p* « 7.SAT in La2ZnIr06, and assuming similar 
ratios of exchanges, Jh/Jk ^ 1, we estimate the dominant 
Jk ~ 15Ar, and Qcw ~ —15AT; however, experiments report 
Qcw ~ —3Ar II 55 I . This discrepancy remains to be resolved. 

In summary, DP Mott insulators are a distinct class of ma¬ 
terials which host strong Kitaev exchange interactions. Our 
study calls for a microscopic understanding of the AFM Ki¬ 
taev exchange, motivates a search for DPs with large T, which 
can stabilize stripes or complex spiral orders. The A-II AFM 
order we find in the AFM Kitaev model is stable against quan¬ 
tum fluctuations for j = 1/2 moments; a detailed study of 
quantum fluctuation effects will be reported elsewhere 1^ . 
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of Physics and Astronomy (SM), and the Scientific User Fa¬ 
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